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(Communicated at the meeting of November 29, 1958) 


Introduction 


In the course of investigations on the role of the imidazolyl group in 
chymotrypsin-catalyzed reactions 1~%) it was deemed necessary to explore 
the possibilities of modifying the imidazolyl group. The reaction with 
iodine seemed to offer a useful way since the iodination can be con- 
veniently carried out under very mild conditions. The literature on the 
reactions of iodine with amino acids and proteins has been reviewed by 
several authors 4-7). 

The reaction of histidine with two moles of iodine in strongly alkaline 
solution leads to the formation of diiodohistidine 8), whereas monoiodo- 
histidine is obtained when one mole of iodine is used 8). The kinetics of 
the reactions of iodine with histidine and with tyrosine in aqueous solutions 
near neutrality have been studied by L1 9). This author found that both 
I, and HOJ, their ratio depending on the pg and the iodine concentration 
used, react with these amino acids. According to Li, the reactions in 
which the amino acids and iodine were present in a mole ratio of 1: 2, 
resulted in the formation of the diiodo derivatives. It was stated that 
the second order kinetics of the reactions were consistent with the reaction 


schemes 


1) D. M. Brouwer, Thesis Leiden, 1957. 

2) D. M. Brouwer, M. J. v. p. Viuer and E. Havinea, Proc. Kon. Ned. Akad. 
Wet. B 60, 275 (1957). 

3) L. Grngaar and D. M. Brouwer, Proc. Kon. Ned. Akad. Wet. B 60, 329 (1957). 

4) R. M. Herriorr, Adv. Protein Chem. 3, 169 (1947). 

5) H. S. Oxcorr and H. FraenKEeL—Conrat, Chem. Revs. 41, 151 (1947). 

8) F. W. Putnam in “The Proteins” (Eds. H. Neurath and K. Bailey, Academic 
Press, Inc., Publishers, New York, 1953) IB, 893. 
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: 2 ; /HOI Fe : 
Tyrosine —¥°!_..- Monoiodotyrosine *“"——> Diiodotyrosine 
and 


Histidine ”291_-+ Monoiodohistidine —"%°-—> Diiodohistidine. 


In these schemes the second steps of the reactions were assumed to be 
much faster than the first ones and, consequently, the monoiodo derivatives 
should occur in very low concentrations only. The validity of this con- 
clusion has been questioned in the case of tyrosine since monoiodotyrosine 
had been isolated from iodinated proteins 4~‘). 

Using radioactive iodine and performing the reactions in strongly 
ammoniacal solutions, Rocue and his co-workers 1°) found that the 
monoiodo derivatives occur as intermediary products in appreciable con- 
centrations. Besides, they observed that, in the presence of an excess of 
iodine, the diiodohistidine formed reacted under these conditions. 

In this connection it should be noted that imidazole is converted into 
tetraiodoimidazole on treatment with an excess of iodine in strongly 
alkaline solutions 1!) whereas on iodination of 2- and 4(5)-methylimidazole 
the triiodo derivatives are obtained). From this, the iodination of 
histidine in the presence of an excess of iodine can be expected to lead 
to the formation of» triiodohistidine («x-amino-f-4(5)-[1, 2, 5(4)-triiodo- 
imidazole ]-propionic acid). In preliminary experiments, in which the 
reactions were followed spectrophotometrically, we found, however, that 
even in solutions near neutrality histidine reacts with more than fou to 
five moles of iodine. Under the same conditions tyrosine reacted in a 
“normal” way, i.e. with two moles of iodine. 


Titrimetric investigation 


In view of the remarkable behaviour of histidine thus found and the 
discrepancies mentioned above, we started to follow the reaction of 
histidine with iodine by means of the same titrimetric technique as used by 
L1%). The mole ratio of amino acid and iodine, however, was 1: 5 *). 

It appeared that at py 6.80 in a 0.1 M phosphate buffer more than 
3.5 out of 5 moles of iodine reacted within six hours; at py 7.40 4 to 4.5 
moles of iodine had reacted within the same time. Qualitatively, the 
dependence of the reaction rate upon py and iodide concentration was in 
agreement with the results reported by Li. Phosphate exerted a strongly 
accelerating influence upon the reaction as found by Li. Borate showed 
a catalytic influence, though in a less degree. Using arbitrary units, the 
following initial velocities were found at py 6.80, the iodide concentration 
being 3.4x 10-2 M: 


10) J. Rocus, 8. Lissrrzky, O. Micuen and 8. Micuez, Ann. Pharm. Franc. 
9, 163 (1951); Biochim. Biophys. Acta 7, 439 (1951). 

Mt) K. Hormann, “‘Imidazole and its derivatives’ (Interscience Publishers Ltd., 
London, 1955) I, 122-124. 

*) For details, see experimental part. 
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0.4 M borate + 0.6 M sodium nitrate: velocity 1 


Ovi +0.7 M a ws 0.9 
0.4 M phosphate +0.6 M a ame Be " 8 
0.2 M = +0.7 M - ee % 4.8 *). 


It was not possible to describe the kinetics of the reaction by a simple 
relation. Assuming that one mole of histidine had reacted with five 
successive moles of iodine, the first step of the reaction being rate de- 
termining, no constant values of k were obtained. With some buffers the 
rate constant showed a considerable drop in the course of reaction time, 
with others it remained constant during the greater part of the reaction. 
When a high phosphate concentration was used, it passed through a 
maximum. It was shown that the decrease of the reaction rate could 
partly be accounted for by a drop of the pH during the reaction that 
occurred even when 0.2 M phosphate buffers were employed. 

The indications, obtained from paper-chromatographic analysis, that 
diiodohistidine reacts with iodine under the present conditions (see 
below), led us to study the kinetics of this reaction. In a 0.1 M phosphate 


ml thio 
40 


3.0 
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lOO 200 300 400 
min 


Fig. 1. Reaction of iodine with diiodohistidine (curve a) and histidine (curve b) 
at pH 6.8 (0.1 M phosphate, 0.7 M nitrate, 0.035 M iodide.) Initial concentrations: 
I, = 4 x 10-3 M, histidine = 0.8 x 10-3 M, diiodohistidine = 0.7 x 107% M. 


*) In separate spectrophotometric experiments with varying concentrations of 
sodium nitrate it was shown that the rate does not depend on the ionic strength 


under the conditions used. 
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buffer, pH 6.80, diiodohistidine appeared to react faster than histidine 
(figure 1). One mole of diiodohistidine reacts with more than two moles 


of iodine. 


Paper chromatographic investigation 


According to Rocns, Lissirzky, MicHEL and MircHen 1°) and Bru- 
NINGS 8), monoiodohistidine is formed in large amounts when the 
iodination is carried out in concentrated ammonia or in a sodium hydroxide 
solution. The possibility should be considered, however, that at considerably 
lower pH monoiodohistidine is present only in very low concentration. 
We therefore examined reaction mixtures at various stages of the reaction 
by means of paper chromatography. The reaction mixtures were prepared 
by treating 0.025 M solutions of histidine with varying amounts of iodine 
which was added in several forms (see experimental part). Analogous 
experiments were carried out with diiodohistidine. 

It appeared that after diiodohistidine had reacted with one mole of 
iodine, its spot could no longer be detected on the chromatograms. A few 
weak new spots and a more intensively coloured one had appeared by 
then which showed up by spraying with ninhydrin as well as with diazotized 
sulphanilic acid. The reaction mixtures obtained with histidine as starting 
material showed a rather weak spot of diiodohistidine as well as a series 
of spots between those of histidine and diiodohistidine. These spots, the 
intensity of which depended on the conditions, appeared almost simul- 
taneously with that of diiodohistidine. Not all these spots corresponded 
to those that were obtained on iodination of diiodohistidine. None of 
the chromatograms showed a spot corresponding to that of monoiodo- 
histidine. 

Interestingly, the R, value of histidine strongly decreased when the 
solution applied to the chromatogram also contained phosphate; other 
anions tested did not have such an influence. The R, value of e.g. 
diiodohistidine was not affected by phosphate. 


Comparison of some imidazole derivatives 


The unexpected behaviour of histidine prompted us to carry out a 
comparative investigation of the iodination of some related compounds. 
Since the preliminary experiments had shown that the kinetic course of 
the iodination could be conveniently followed by the decrease of the 
optical density in the visible region, the further measurements were 
performed in this way. This method offered the additional advantages 
that, because of the lower concentrations of the reactants, the pH remained 
more constant and that the concentration of iodide, which decreases the 
reaction rates, could be kept lower. The concentrations generally used 
were about 8.0 10-4 M for iodine, 1.0 10-4 M for the other reactant 
and 2.4 10-3 M for potassium iodide. 
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Figure 2 presents the results of a few experiments in a 0.092 M phosphate 
buffer of pH 7.40. It is seen that one mole of imidazole reacts with four 
moles of iodine and one mole of 2-methylimidazole with three moles of 
iodine. In both cases a precipitate formed. In view of the insolubility of 
1, 2, 4, 5-tetraiodoimidazole and 1. 4. 5-triiodo-2-methylimidazole it seems 
safe to conclude that under the present conditions the same products are 
formed as in strongly alkaline solution (cf. 1.c.12)), 
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Fig. 2. Reaction of iodine with histidine (curve a), 2-methylimidazole (curve c) 

and imidazole (curve b) at pH 7.4 (0.092 M phosphate, 0.0024 M iodide). Initial 

concentrations were 8.0 x 10-4 M for iodine and 1.0 x 10-4 M for the other 
compounds, respectively. 


Histidine reacts rapidly with five to six moles of iodine, whereupon 
iodine continues to be consumed in a very slow reaction. A definite end 
point was not reached. 

Whereas in the case of histidine the reaction is strongly accelerated by 
phosphate — the initial rate being about fifty per cent lower in a 0.046 M 
phosphate buffer of the same pH (7.40) and comparable ionic strength — , 
the rates of iodination of imidazole and of 2-methylimidazole were found 
to be independent of the phosphate concentration. 

Evidently, neither the drastic course of the iodination of histidine nor 
the catalytic effect of phosphate are features characteristic of imidazolyl 
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compounds in general. In order to check the alanyl part of the molecule 
as to its possible reaction with iodine, the behaviour of alanine and 
tyrosine under the conditions used was investigated. Alanine did not 
react, tyrosine reacted with exactly two moles of iodine, the reaction, at 
pH 6.00 and at pH6.75, being much faster than that of histidine under 
the same conditions. The latter results are in agreement with the data 
reported by Li and by Rocue et al. 

The iodination of diiodohistidine is also accelerated by phosphate. In 
0.092 M phosphate of pH 6.74 the reaction proceeds about 20 % faster 
than in 0.046 M phosphate/0.046 M sodium nitrate. More than four moles 
of iodine had reacted within three hours. Histidinemethylester and 
histamine react in essentially the same way as histidine as can be seen 


from figures 2 and 3. 


100 200 1000 2000 3000 
min 


Fig. 3. Reaction of iodine with a) histamine, initial concentration 1.0 x 10-4 M, 
at pH 7.35 (0.046 M phosphate, 0.046 M nitrate). 6) benzoylhistidine, initial 
concentration 0.9 x 10-4 M, at pH 7.5 (0.092 M phosphate). c) acetylhistamine, 
initial concentration 1.0 x 10-4 M, at pH 7.40 (0.092 M phosphate). Initial iodine 
concentration 8.0 x 10-4 M and iodide concentration 2.4 x 10-3 M in all experiments. 


The iodination of «-N-benzoylhistidine (and of «-N-acetylhistidine), 
however, seems to follow a different kinetic pattern (figure 3), These 
compounds reacted readily with two moles of iodine, the first mole of 
iodine being consumed much more rapidly; then a very slow reaction 
followed. Indications were found of a catalytic influence of phosphate. 
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No coloured products were formed during the reaction in contradistinction 
to what happened with diiodohistidine (cf. experimental details). When 
the unchanged iodine was removed by addition of some tyrosine, the 
solution became colourless. 

a-N-acetylhistamine behaved in an analogous way (figure 3). After two 
moles of iodine had reacted rapidly, the reaction stopped. Phosphate 
catalyzed the reaction. 

N-benzoylglycine was found not to react with iodine. 


Discussion 


From our experiments and from the literature data cited it follows 
that the formation of monoiodohistidine in appreciable concentrations 
depends largely on the conditions under which the iodination is performed. 
Evidently, the rates of substitution by the first and the second iodine 
atom are differently influenced by the reaction conditions, the first reaction 
being more strongly accelerated by high alkalinity of the medium. It seems 
attractive to explain the different results reported in the literature 4-7, 9. 10) 
with respect to the occurrence of monoiodotyrosine as a product of the 
iodination of tyrosine on a similar basis. 

From the chromatographic investigations it follows that some of the 
products resulting from iodination of histidine, are not formed as over- 
iodination products of diiodohistidine. Besides, it depends on the conditions 
used to which extent these products are formed. 

As stated before, the behaviour of histidine towards iodine is neither 
inherent in the imidazolyl group nor in the «-amino acid grouping as such. 
Apparently, some cooperation of the imidazolyl group and another part of 
the molecule (intramolecular catalysis?) is essential to the drastic course 
of the reaction with iodine. Since histamine and histidinemethylester react 
in the same way as histidine, the conclusion seems justified that the 
carboxyl group is not important in this connection (although there is a 
slight difference between the course of reaction of «-N-acetylhistamine 
and that of «-N-acetylhistidine or benzoylhistidine (Comp. Fig. 3)). 
The marked change in the reaction course that occurs on acylation of the 
amino group — both with histidine and with histamine — suggests that this 
group has to do with the remarkable behaviour of these compounds. 
The catalytic influence exerted by phosphate probably should be 
attributed to an interaction between the catalyzing ion and histidine, 
histamine, tyrosine, etc., because such an influence is absent in the cases 
of imidazole and 2-methylimidazole. A direct relation between this 
catalytic effect and the remarkable course of the iodination of histidine 
and histamine seems questionable since the catalytic effect of phosphate 
is also found with «-N-acylhistidines and with tyrosine °). 
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Titrimetric measurements 


50 ml of a 1.6 x 10-3 M solution of histidine or diiodohistidine, respectively, 
and an equal volume of an 8.10-* M solution of iodine were mixed in a glass-stoppered 
flask. The solutions were freshly prepared in the buffer chosen *) and thermostated 
at 25° before mixing; the iodine solution contained also the desired amount of 
potassium iodide. At predetermined intervals, 5.0 ml aliquots were acidified with 
10 ml of N hydrochloric acid and the unreacted iodine was titrated with 0.01 N thio» 
starch being used as indicator. In the case of diiodohistidine a yellowish reaction 
product formed after some time. Therefore, the indicator solution was added when 
slightly less than the expected equivalent of thio was added. 


Paper chromatography 


0.025 M solutions of histidine and of diiodohistidine were treated with iodine 
that was either dissolved in a potassium iodide solution, in ether or in ethanol, 
or added as KI + H,O,. The pH was kept constant by a buffer or by continuous 
addition of base. 

10-20 wl samples of the reaction mixtures were applied to Whatman no. | filter 
paper after the reaction had been stopped by addition of acid or by extraction of 
the excess of iodine with ether. The chromatograms were developed by the ascending 
technique for sixteen hours with butanol-acetic acid-water (78: 5:17 or 40: 10: 50) 
as the eluent. Duplicate chromatograms were sprayed with ninhydrin solution 
prepared according to BorssonNaAS !”) and with diazotized sulphanilie acid, respec- 
tively. The latter reagent was conveniently prepared by mixing, in the given order, 
4 ml of a 0.5 % solution of sulphanilic acid in 1 % hydrochloric acid, 1 ml of a 4 % 
solution of potassium nitrite and 5 ml of a 10 % sodium carbonate solution. One 
ug of histidine was easily detected with the aid of this reagent. 


Spectrophotometric method 


2.0 ml of a solution of iodine containing potassium iodide was added to 25.0 ml of 
the buffered solution *) of the product to be iodinated. The final concentrations usually 
were 8.0 x 10-4 M iodine, 2.4 x 10-8 M potassium iodide and 1.0 x 10-4 M product. 

An aliquot was immediately transferred to a 10 mm glass-stoppered cell and the 
optical density was measured as a function of time at 450 my in a Unicam SP 600 
spectrophotometer. The measurements were carried out at room temperature. 

In the cases of imidazole and of 2-methylimidazole a precipitate (see text) 
settled on the walls of the reaction vessel and clear samples could be taken after 
some time. 


Summary 


Treatment of imidazole and of 2-methylimidazole with iodine in neutral 
aqueous solutions results in the substitution of the nuclear hydrogen 
atoms. One mole of histidine, however, reacts rapidly with as much as 
five or six moles of iodine, which process is followed by a slow further 
consumption of iodine; histidinemethylester and histamine react in a 


*) As a rule, the buffer solutions also contained sodium nitrate so that the 


buffer concentration could be varied without considerable change of the ionic 
strength. Compare also the footnote of page 95. 
2) R. A. Borssonnas, Helv. Chim. Acta 33, 1975 (1950). 
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similar way. Monoiodohistidine is not formed in detectable concentrations. 
In the particular case we investigated the reaction of iodine with diiodo- 
histidine proceeded more rapidly than that with histidine. Not all of the 
iodination products of histidine are formed via diiodohistidine. Under the 
conditions used acetylhistamine and benzoylhistidine react with but two 
moles of iodine (a slow reaction seems to follow in the case of benzoyl- 
histidine) ; tyrosine reacts with two moles of iodine; alanine and benzoyl- 
glycine do no react at all. 

Phosphate accelerates the iodination of histidine, histamine (and 
benzoylhistidine) but not that of imidazole and 2-methylimidazole. 

A discussion of the phenomena is given. 


Laboratory of Organic Chemistry, 
The University, Leiden 


PHYSICAL CHEMISTRY 


PRECISION VISCOMETRY OF POLYISOBUTYLENE IN TOLUENE 
BY 


H. W. ZWEEGMAN anv J. J. HERMANS 


(Communicated at the meeting of November 29, 1958) 


The viscosity of solutions of polyisobutylene in toluene, ranging in concentration 
from 1.2 x 10-5 to 1.7 « 10-3 g/ml is found to be Newtonian for shearing stresses 
below about 5 dyne/em?. There are indications for the occurrence of adsorption, but 
it appears that the layer of adsorbed polymer cannot be removed by rinsing with 
the solvent. This may be an explanation for the fact that the relative viscosities 
derived from the efflux times are independent of the capillary diameter. The curve 
of 7,,/¢ versus c shows pronounced downward curvature below a concentration 
of about 10-4 g/ml. 


1. Introduction 


The viscous behavior of very dilute polymer solutions has been the 
subject of numerous recent studies. According to STREETER and BoYER 
(1951), CurLeR and Krwpaui (1951), Moore and RussEuut (1952) and 
GIESEKUS (1954), the curve of 7s)/c versus c shows pronounced downward 
curvature when c is decreased. In later work, STREETER and BoyEr (1954), 
BaTzeR (1954), Kapapra (1954, 1955), and Takepa and Enpo (1956) 
reported a minimum in the curve, sometimes followed by a maximum. 
No such anomalies were found, however, in the work of RoTHMAN, SIMHA 
and WEISSBERG (1950) and of Kooy and Hermans (1955). More recently, 
TUYNMAN and HERMANS (1957) found anomalies which could be explained 
partly by adsorption and partly by non-Newtonian effects, but the non- 
Newtonian character of the flow has not been confirmed by the recent 
work of OuRN (1958a). 

OuRN (1954) was the first to discuss the importance of the adsorption 
of the polymer which, on the one hand, results in loss of polymer and, 
on the other hand, decreaes the effective capillary diameter of the visco- 
meter. A discussion of the anomalies reported in terms of the so-called 
“critical concentration,” i.e. the concentration at which the solute mole- 
cules begin to lose contact with each other, can be found in several of the 
articles quoted; we mention, in addition, Kawatr and Sarro (1957), Hrrat 
(1957) and KriaBaum (1957). 

No attempt has been made in this introduction to give an exhaustive 
review of the literature, because this has been done admirably by Ourn 
(1958a), who was able to show that almost all the anomalies disappear 
when appropriate corrections are made for adsorption effects. Compare 
also the discussion of Umstiitter’s paper (1958) by Ourn (1958b). It is 
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only at extremely low concentrations, where the results become less 
reliable, that certain indications for a downward bend of Nop[C Versus C 
are found. For polystyrene of molecular weight 8105 in toluene, this 
occurs at c about 10-5 g/ml. The results reported below indicate a similar 
downward curvature for polyisobutylene, M— 1.6 Xx 106, in toluene below 
a concentration of about 10-4 g/ml. 


2. Hxperiments 


The polyisobutylene was a sample of Vistanex L 300, kindly put at 
our disposal by Goodrich Research Center, Brecksville, Ohio. It was 
dissolved in pure cyclohexane and filtered through a P3 glass filter to 
remove gel. The solution was sprayed through an atomizer into a mixture 
of 4 1 methanol and 2 | acetone, stirring vigorously. The precipitated 
polymer was dried and then ground to a powder at liquid air temperature. 
Light-scattering studies in cyclohexane, using a Zimm plot, revealed that 
the weight average molecular weight was 1.6106, while the radius of 
gyration was 870 A. For the effect of the polymer on the refractive index 
n we found, in cyclohexane, dn/de=0.10 ml/g. The intrinsic viscosity in 
toluene at 25°C was 325 ml/g. 

The toluene (Merck, p.a.) was made dustfree by repeated filtration 
through P3 glass filters. All glassware was made dustfree by Thurmond’s 
method (1952). The solutions were prepared by repeated dilution of a 
stock solution, and the loss of polymer by adsorption was minimized by 
the procedure described by TuynmMAan and HERMANS (1957). 

The temperature control and the technique of measuring the efflux 
time have been described in previous papers by Kooy and HERMANS (1955) 
and TuyNMAN and HeRMANS (1957). The viscometer (Fig. 1) was of the 
suspended level type (Ubbelohde) and contained three spherical reservoirs, 
each with its own capillary. The dimensions are given in Table I. 


TABLE I 


Dimensions of the viscometer, and the values of the constants A and B of Eq. 1 

for toluene at 25° C. Capillary radius R and capillary length L in cm, bulb volume 

V in em3, A in dyne/em?, B in dyne sec/em*. The standard deviations o4, and o, 
are likewise in dyne/em? and dyne sec/em? respectively 


renee ee re err eee ee 


Capillary R L V A O4 NOH 218; 10-56, 
ee _ 
1 0.0318 110.5 0.417 4326 9 6.419 0.011 
2 0.0219 68.0 0.300 4489 9 12.700 0.024 
3 0.0188 45.0 0.184 4355 10 9.377 0.021 


The pressure which was applied to study the effect of the shear-rate 
was measured in a manometer which allowed readings to within 0.01 mm, 
corresponding to better than 1 dyne/cm?. This manometer was perfectly 
reproducible and showed no hysteresis. Fluctuations in the pressure were 
very small. Yet it appears, from the experimental results, that the pressure 
measurements were not quite as satisfactory as they seemed to be (see below). 


LO5 


In all our experiments the Hagenbach correction for kinetic energy 
was less than 0.07 %. Since this correction is made both for the (very 
dilute) solution and for the solvent, its effect on the relative viscosity is 
extremely small so that it is not important to know the exact value of 
the Hagenbach constant. 

The effect of the change in hydrostatic head during the measurement 
has been calculated by us for a spherical reservoir. For pure toluene it 
amounts to 0.16 °% when there is no external pressure difference, and it 
decreases rapidly with increasing pressure. Here again, the effect on the 
relative viscosity is, of course, much less than on the viscosity itself. 
Finally, the drainage error can be shown to be negligible. 

Our procedure consisted of first correcting the pressure for the effects 
of kinetic energy and change in hydrostatic head. Then, if A is the effective 
hydrostatic head in the absence of an external pressure, the relation 
between the corrected external pressure, »,, and the efflux time ¢ for a 
Newtonian liquid of viscosity 7 is 


(1) pet A=Bit 


where b= 87LV/zR*. The value of A for toluene was determined statisti- 
cally, using 1/f as independent and p, as dependent variable. The result 
is given in Table I, where o4 and og represent the standard deviations of 
A and B. 


3. Results and Discussion 
The total pressure difference p and the maximum shear stress t follow 
from the relations 


(2) P=pet+A; t=pR/2L. 


Plots were made of pt versus p for six solutions with concentrations 
ranging from 1.70 10-3 to 1.26 10-5 g/ml. These plots did not deviate 
systematically from a horizontal line, which means that in the range of 
shear stresses investigated the solutions were Newtonian. There was, 
however, an unexpected scatter about the horizontal line, and this scatter 
was systematic in the sense that for each of the solutions, the deviation 
for a given pressure had the same sign and approximately the same 
magnitude as that found for toluene at the same pressure. ‘This is 
demonstrated by the example given in Fig. 2: the open circles indicate 
the values of pt/Bo where Bo is the average value of pt for pure toluene; 
the closed squares are the result of dividing each pt for the solution by 
the pt-value for toluene at the same pressure. We believe that the results 
shown can be explained only by assuming that there was some unknown 
error in our pressure measurements. Part of the difference in pt/(pt)o 
between the circle and the square at the lowest p-value, where pe=0, 
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Fig. 2. Relative viscosity versus pressure for a solution containing 0.033 x 10-% 
g/ml polyisobutylene in toluene at 25° C; data obtained in capillary nr. 3. 


could perhaps be explained by a failure to account accurately for the 
change in hydrostatic head during the measurement, but this explanation 
cannot possibly apply to the other points. 

The results obtained are such that the error in p does not affect our 
conclusion that 7 is independent of p; nor does it affect the value of the 
relative viscosity 7, to any appreciable extent. It was found, further, 
that 7, was the same within experimental error for all three capillaries. 
Yet, we do have an indication for adsorption of the polymer by the wall 
of the capillary: the efflux time for toluene was slightly lower before than 
after measuring that for the solutions, even though the viscometer was 
rinsed repeatedly with toluene after being filled with the solutions. It 
would seem, therefore, that the rinsing did not remove the adsorbed 
polymer, perhaps because toluene is a comparatively poor solvent for 
polyisobutylene. This view is confirmed by experiments by FENDLER et al. 
(1956). From the difference in efflux times for toluene before and after 
the solution series, we calculated a thickness of about 500 A for the 
adsorbed layer, but this is no more than a rough estimate. The solutions 
were examined in the order of descending concentration, starting with 
c=1.70x 10-3 g/ml. The fact that we found values of 7s, which were 
independent of the radius of the capillary used could be explained by 
assuming that the adsorbed layer was formed when the viscometer was 
filled with the most concentrated solution and remained essentially 
unchanged throughout the series of dilutions. 

Fig. 3 shows that the experimental values of 7,,/c are anomalous below 
about 10~4 g/ml. This effect is either real or else due to sources of error 
which we have failed to detect. It cannot be due to non-Newtonian 
behavior, or to the effect of adsorption on capillary diameter, because the 
measurements themselves indicate that such effects were non-existent. 

Finally, since Fox, Fox and Fiory (1951) have found that polyiso- 
butylene in various solvents shows non-Newtonian viscosity at shear-rates 
higher than those examined by us, a series of experiments were carried out 
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Fig. 3. Reduced specific viscosity versus concentration of polymer. 


at shear stresses ranging from 8 to 22 dyne/cm2. These confirmed the non- 
Newtonian behavior at shear stresses larger than 8 dyne/cm2?. Below this 
value the curve of viscosity 7 versus shear stress t becomes horizontal, in 
accordance with the requirement that 7 must be an even function of rt. 
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PHYSICAL GEODESY 


THE OUTSIDE GRAVITY FIELD UP TO GREAT DISTANCE 
FROM THE EARTH 


BY 


F. A. VENING MEINESZ 


(Communicated at the meeting of December 20, 1958) 


Flying missiles and Sputniks have suddenly increased the importance 
of knowing the gravity field of the earth up to great distance. Although 
the writer formerly published formulas about this matter 1), it may at 
present be useful to go into more detail. 

We start from the formula for the outside potential of a spheroidal 
earth; the coordinates used are the radius 7 to the earth’s center of 
gravity and the geocentric latitude g. We shall assume that the outside 
point where we investigate the potential is not rotating round the earth’s 
axis. Up to the fourth order term, i.e. up to the second order of the earth’s 
flattening, we have 


(1) U=ke [1 = al (5 sin? y-3)+z (sins pz sin? +35) val 
where k is Newton’s constant, J the mass of the earth, C the moment of 
inertia round the rotation-axis, H the “mechanical ellipticity’’, i.e. the 
ratio (C—A)/C, in which A is the moment of inertia round an axis in 
the equator plane, and finally, D a higher order moment of inertia 
divided by M, given by 


M 
1 1225 3 Grrr 3 
(2) D=z{ 7 ae (sint y—= sin p+) dm. 
0 


We know that, because of its rotation, the earth in its general shape 
assumed the figure of a flattened spheroid, and research by Clairaut, 
George Darwin, Willem de Sitter and others have shown that the equi- 
librium figure of a rotating fluid body of the same volume, mass and 
general density distribution as the earth is a spheroid closely approaching 
a rotation ellipsoid; both surfaces do not deviate more from each other 
than about two meters. It is clear, therefore, that the geodesists, in taking 
a flattened ellipsoid as reference figure for the earth, are near to the true 
figure. 

Examining now our problem in greater detail, and taking for the 
figure of the earth the geoid, which is the equipotential surface of the 


1) HEISKANEN and VENING MEINESZ, “The Earth and its Gravity Field’, 
3-3, 3-4, 3-5, 3-6, 3-9, McGraw-Hill, New York, 1958. 
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earth (therefore everywhere at right angles to gravity) at mean sea-level, 
we find that the real earth in three ways deviates from the rotation- 
ellipsoid; all three, however, only slightly affect the gravity field at 
greater distance from the earth. In the first place the equilibrium spheroid 
is not exactly the rotation-ellipsoid. In the second place the geoid is not 
precisely coinciding with the equilibrium spheroid. In the third place 
the real surface of the earth’s masses shows topography and, therefore, 
does not coincide with the geoid. In these masses we obviously have 
to include the water masses at the earth’s surface; we can, however, 
show that the masses of the atmosphere are too small to affect our 
problem appreciably. 

Of the three deviations mentioned the first is evidently likewise negli- 
gible. For the third this is also true; as the rigid crust of the earth, of a 
thickness of about 35 km, is floating on the heavier and plastic subcrustal 
layer, the masses of the topography are, therefore, compensated by 
masses of opposite sign lower down. This phenomenon, called isostasy, 
makes for the total mass in every vertical column of the same horizontal 
cross-section and reckoned up to the same depth being about the same 
over the whole earth’s surface. The differences in location of the masses 
in each column are too small for appreciably affecting the outside gravity 
field at greater elevation than some fifty kilometers. Near the continental 
coasts, which constitute the largest topographic features at the earth’s 
surface, we have perhaps to increase this limit up to a height of 100-200 
km, but the deviations in figure of the equipotential surfaces and in 
gravity are local and average out over a belt along the coasts of a few 
hundred kilometers broad. 

The largest deviations are caused by the second of the three causes 
mentioned above. The masses of the earth appear to deviate sufficiently 
from the regular layering which would correspond to fluid equilibrium 
for causing the geoid to deviate from the equilibrium spheroid over 
horizontal distances up to one thousand kilometers or more and over 
vertical distances up to 40 meters. Similar deviations must be present 
in outside equipotential surfaces and we shall have to study the vertical 
distances above the earth’s surface at which they become negligible. 
We shall take this up at the end of this paper. Compared to the dimensions 
of the earth as a whole they are of limited extent. 

Replacing the equilibrium spheroid by a rotation ellipsoid of an equa- 
torial radius a and a flattening «, we can derive that for this purpose 
we have to put 


(3) nao (ga-se/) 


where c’ is a coefficient dependent on the velocity w of rotation of the 
earth in the way given by 


4 1 _ wras 
(4) =. 


1a 


It can likewise be shown that the coefficient before the brackets of the 
second term between the square brackets of formula (1), which as a factor 
contains the ‘mechanical ellipticity” H, is a function of the flattening « 
and the coefficient c’ in the way given by 


2a —c' — 202+ 2xc’ +e5 = 
Le 


( CH CH _@ 5 (2 sey a2 

(5) , e , : 
re oem re ra ac ee 

‘ 3 (20 C oP +5 ae!) &. 


Our formula (1) now gives for the outside potential of an ellipsoidal 
earth at a point at a distance r from the earth’s center and a geocentral 
latitude 


M 1 i j 9 Ors 1 
= 7.74 Ae (Oey eee ea Dae cena 
ve k= |1 5 (20 C x? + = ae’)  (Fsin® g 5) + 


7 i) 
+ (5 02 5 ac’) S (sinty —S sin? g + 2 a) en | 


(6) 


For making this formula complete for the rotating earth, we have to 
add the potential of the centrifugal force, which is given by 


(7) U, = 4w? r2 cos? 
In this formula w is the velocity of rotation; it can easily be computed. 


The radius r, of an arbitrary point of the earth-ellipsoid, which has a 
great axis 2a and a flattening «, is up to the second order of « given by 


(8) re=a[1—(x+5 a2) sin? p +5 a%sint gy... ]. 


Because the earth-ellipsoid is supposed to be an equipotential surface of 
the rotating earth, this value of r,, introduced in the sum U + U, must 
give a result that is independent on the geocentric latitude gy. This condition 
provides us with the equations (3) and (5) already mentioned. 

We shall now use our formula (6) for the earth’s gravity field for points 
not rotating with the earth. The formula for the radius of equipotential 
surfaces of this field must up to the second order of their flattening have 
the character of the radius of a spheroid. We shall indicate this radius 
by o and the equatorial radius of an external equipotential surface by 
a,; we shall give the same subindex to the flattening and to the corre- 
sponding 4th order term. We write, therefore, for the radius of the external 


equipotential spheroid 
(9) @ = a.[1 — (a, + %4) Sin? V+ oy sin* ¢] 
Introducing this value of r in equation (6) we obtain two equations 


expressing that the value of U must be independent of ¢; the coefficients 
of sin? and sin’ must, therefore vanish. This give us 


1 Seat ae! a ae! a) # 
(10a) oy = 5 (24-0! a8 + 5 ae) B35 (—a Saad mmc. ig 
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(10b) a= 5 (Ba — axe’ —c'?) 
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From formula (10a) we can draw the conclusion that the first term of 
the flattening «,, which forms by far its main part, is inversely proportional 
to the square of the distance a, from the earth’s center to the equator 
of the equipotential surface. We could also have drawn this conclusion 
from the fact that the flattening constitutes a second order spherical 
harmonic of the earth’s gravity potential; divided by gravity (inversely 
proportional to a2) and by a,; the second order spherical harmonic of 
the gravity potential itself must be inversely proportional to ap. 

For « and c’ we may now introduce the values for the internationally 
adopted earth-ellipsoid: 


(11) x = 0,0033670 = 1/297.00 —c’ = 0,0034614 
We thus obtain 
(12) % _ 0.48654 % + 0.00035 © 

“ % a" 


which shows that the 4th order term is practically negligible. We, there- 
fore, can state that the flattening is inversely proportional to the square 
of the equator-radius of the equipotential spheroid. We get the following 
values 


a 1 2 3 4 
(13) 
[o9 
* 


£= 0.48688 0.12166 0.05406 0.00487 


We see that the suppression of the rotation makes the equipotential 
spheroid with the same equator as the earth-ellipsoid have about half 
of the flattening of the earth-ellipsoid. We also see the quick decrease 
of the flattening with the distance. 

We shall now turn to the outside gravity field itself. Differentiation 
of U (formula 6) with respect to 7 gives us the gravity component g, in 
the direction of this radius and with respect to m the much smaller 
gravity component gs in the south direction at right angles to the radius. 
We obtain 


oU M , 9 a’ /3 . 1 
SF =a (eet fae) $ Gat) 
us \ a < 2 Sl — 7 HC ) (5 Sin’ p—5) + 
7 5 W\ es. : 
+5 (50% — 3 ac’) S (sint p— 2 sin? g + =) 
oU M ae os 
(14b) 2 ae =z (2a —¢') 5 sin y cos g. 


As, obviously, the effect of the south component gs on the total value of 
gravity is of the order of the flattening, we have not further derived the 
formula for gs than up to the first order of «. Under these circumstances 
we can determine the resultant of both components by means of the 
formula 


9 oft Q] 
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and obtain for g an accuracy up to the second order of the flattening ; 
by replacing sin? y cos? p by sin? y—sin4 y we get 


M ; 5 NCE Bc 
(9=h | 1—(20—c —a2+ Fac") (5 sin? yp —3) + 
(16) i Ls 
7 SINCE os Om om, 3 1 Asey oben ye ot : 
+5 (50°30!) © (sint p—5 sin? p+ 35)+5 (20 —c')? & (sin? p—sin‘ g) | 


By dividing gs by g, we can obtain the tangent of the angle 6 between 
the direction of gravity and the geocentric radius r and because this angle 
is of the order of the flattening. we can replace the tangent by the angle 
itself. By multiplying the ratio by @” we obtain this angle in seconds of 
arc. We find 


een lias c , a? A B ; az 3 : l 

(17) @ =0" (2a4—c’) =z SID P Cos p E + (206!) = 5 sin? y—5) |. 
Introducing the values of « and c’ of the earth, the formulas (16) and 
(17) become 

M az/3 . il 
g=k= 1 — 0.972 «= (5sin?p—s)+ 
(16A) \ r (5 3) 
ae 650 a2 5 (sin! —§ sine +5) +0 472 «2 = (sin? » —sin4 ) 
\ T “re rt Y 7 YP 35 . a r (s Ost 7 | 


" ” a. Ch? es if 
(17A) 6” = 198420 x & sin pcos g [1+0.972 « 5 (5 sin? y—5)]. 


For obtaining the direction of gravity the angle 6 has to be added 
to g. We thus have derived gravity outside the earth in magnitude and 
direction. 

We shall now take up the problem how we can determine the relatively 
slight deviations in the equipotential surfaces and in gravity caused by 
the disturbing masses which bring about the gravity anomalies and the 
deviations between earth-ellipsoid and geoid. Their effect is limited 
compared to what has hitherto been considered. It is, therefore, only 
appreciable up to small distances outside the geoid. 

We apply Stokes’s method, adapted to this case, and consider the 
undisturbed equipotential surfaces at geoid level and at an arbitrary 
distance outside it as spheres with radii a and a,. We shall denote the 
gravity anomalies on these surfaces by 6g and 6g, and the distances 
between the disturbed and the undisturbed surfaces by N and N,. 
Denoting a/a, by u we have w<1. By means of Stokes’s method we can 


write 
S=470a* 


(18a) N=, [ F,égds 
0 


_ 479.0, 


with 


(18b) F, =2 +——5 cos y—32—3 cos yn } (1—u cos p + ») 


where 


(18c) v = (l—2u cos p+ u?)* 


These formulas give us the disturbance in the outside equipotential 
surfaces caused by the disturbing masses in the earth which cause the 
gravity anomalies. 
For deriving the gravity anomalies on these outside equipotential 
surfaces, we can write 
S=4na' 


1 = Y 
= 4 (2 ta unre ds. 
(19) bg ina J q [ a n+1)Pn ] 


In the following way we can derive the quantity between the square 
brackets and obtain the solution. 


G, = ne (2n+1) Py unt? = 2u%t 1S (nt 4) Phu * = 2uth 5 (> Pawn), 
This gives 
d 
(20a) G= 2083, (eo) 
where 
(20b) q= > Pnau® = (1—2u cos p+ u?)-*—1— cos yp. 
2 


We thus obtain the solution (19) in the shape 


S=42a* 


si [ G.ogas. 


(20c) 69, = rer 


0 


Although we thus can derive the outside equipotential surfaces and 
the gravity field, this solution does not provide us with a quick insight 
in the distance from the earth’s surface up to which the earth’s disturbing 
masses have an appreciable effect. We can better solve this problem by 
making use of the fact that we can represent the outside gravity anomalies 
corresponding to anomalies on the geoid distributed according to a 
spherical harmonic Y, of order n by the formula Y,wu"+2. We thus 
can derive that the gravity anomalies sink to about one percent at an 
elevation of 60.000/n km. For the belts of strong negative anomalies in 
island-are areas, as e.g. in the areas of the Indonesian and Caribbean 
archipelagos, this is the case for elevations of 400-500 km. This is a good 
illustration of the relatively small distances from the earth to which the 
effect of disturbing masses in the earth is limited. 


GEOPHYSICS 


THE RESULTS OF THE DEVELOPMENT OF THE EARTH’S 
TOPOGRAPHY IN SPHERICAL HARMONICS UP TO THE 
3lst ORDER; PROVISIONAL CONCLUSIONS 


BY 


F. A. VENING MEINESZ 


(Communicated at the meeting of January 31, 1959) 


The Netherlands Organization for Pure Research has made possible 
an important undertaking, the development in spherical harmonics up 
to the 3lst order of the earth’s topographic elevation above sea-level. 
It was done by the Mathematical Centre in Amsterdam, where Dr. Ir. D. J. 
Hofsommer and Mr. M. L. Potters carried out the scientific research and 
supervised the computations. This study and work will be the subject of an 
extensive publication. Prof. Ir. G. J. Brutys of the Geodetic Laboratory 
of the Polytechnical University at Delft with his collaborators got together 
the 40680 figures of the topographic elevation — about one for each square 
degree —required for this great work. The writer feels deeply indebted 
to all who collaborated for this important achievement. 

The whole undertaking was a repetition and extension of a similar 
development in 1919-1922 by Prof. Dr. A. Prey of the University at 
Prague; he carried it out up to the 16th order. The 289 coefficients which 
Prey obtained did not appear to show much regularity or system, but 
for a reason which he shall presently mention the writer in 1950 combined 
all coefficients of the same order in one figure, representative for that 
order, and thus found 16 figures forming a curiously regular curve. It 
implied interesting conclusions 1). 

There was reason to repeat the development and to extend it to a 
higher order, because since 1922 much more became known about the 
topography of the ocean-floors; the echo-sounding method provided an 
enormous number of data. Secondly an extension was desirable because 
the conclusions drawn from the first development promised further 
important results if it could be carried out to greater detail. 

The middle curve of fig. 1 shows the results of the new development 
which we shall denote by t. Each ordinate of this curve gives a repres- 
entative figure for the order n as indicated on the horizontal axis. As 
representative figure was chosen the root of the mean square, taken 
over the whole earth’s surface, of the topographic elevation represented 


1) A remarkable feature of the earth’s topography, origin of continents and 
oceans, Proc. Kon. Ned. Akad. v. Wetensch., Amsterdam, Ser. B, 14, 3 (1951), 
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by the 2n+1 sub-terms of the term of the nth order. The formula for 
this nth order term, as well as that for the mean square of the elevation 
represented by it, are given in formulas (2) and (3) on page 4 of the 
writer’s publication mentioned in the footnote. The results ¢ of the new 
development are also given by column 3 of Table I. 


TABLE I (J, t, s in meters) 


n l t 8 DL fh S n l t Ss L gp S 
a 
1 1205 979 869 1704 1385 1229] 16 326 239 192 5376 3942 3166 
2 1065 841 729 2610 2060 1786 1 47 275 239 6074 4810 4181 
oF L148" 905 7833979 8135 27 54 212 192 4699 3921 3551 
4 


1141 875 738 5104 3913 33C 
5 1109 815 657 6072 4464 3599 | 20 


6 694 509 409 4500 3299 2651 | 21 249 192 163 5346 4127 3504 
7 541 481 454 4049 3600 3398 | 22 214 175 156 4818 3937 3509 
8 416 338 300 3532 2868 2546 | 23 202 143 110 4741 3360 2584 
9 488 384 332 4627 3643 3150 24 191 163 150 4676 3993 3674 
10 484 373 316 5079 3912 3314] 25 184 156 143 4684 3977 3646 


11 378 289 243 4344 3320 2792 | 26 244 186 156 6465 4928 4133 
12 316 247 212 3940 3085 2648 | 27 182 135 110 5001 3712 3024 
13 343 271 235 4627 3656 3170] 28 168 135 119 4776 3847 3391 
14 341 293 271 4936 4246 3926 | 29 161 127 110 4749 3746 3245 
15 305 243 212 4731 3765 3284 | 30 210 156 127 6413 4757 3873 


31 101 101 4101 3181 3181 3181 


5 
84 156 143 3583 3041 2788 
33 186 163 4767 3812 3340 


Besides this main development, a second one was made ~—as also Prey 
had done—of all the topographic elevations of negative sign, putting 
all those of positive sign zero. The results for the 31 representative figures 
thus obtained, which we shall denote by s, are given by the lower curve 
of fig. 1; they are also listed in column 4 of Table I. This second develop- 
ment, which we may indicate as the ocean-floor development, makes it 
possible to use our results for solving problems in which the ocean-depth 
is not simply used as negative elevation in the same way as continental 
elevation is used as positive elevation. As an instance of such a problem 
we may mention the computation of the thickness of the chemical crust, 
when basing ourselves on the supposition of isostatic equilibrium. 

From the two lists of representative results the writer derived the 
representative figures for each order n of a development in spherical 
harmonics of the positive topographic elevations, putting zero all the 
negative values. The representative figures of this development — which 
we may indicate as the development of the continental topography, and 
which we shall denote by /—are given by the upper curve of fig. 1 as 
well as by column 2 of Table I. These values were derived from the 
values of t and s by means of the following formula, in which ¢ is the 
ratio of the total surface of the globe to the land surface. 


(1) 1 = Vt#2+(c—1)(#2—8?). 
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As the land surface is 29.17 % of the whole earth’s surface, we have 
c= 3.428. 

Studying the values of ¢ and s we see that up to the 16th order they 
check well with the values of Prey, though a few deviations may be 
noted. The most striking ones are larger values for »=7 and smaller 
ones for n= 15. The general character of diminishing values with increasing 
n remains true, also for the continuation from n=17 to n=31. For the 
further study of the development it is, therefore, useful to multiply all 
the values, in the same way as it was done for the publication of 1951, 
by n, or better still by nt(n+1)* as, for theoretical reasons, it was 
considered preferable. Doing so for the three quantities 1, t, and s we 
obtained the values L, 7’, and S given by the last three columns of the 
table and represented in fig. 2. We see that in great lines the three curves 
are horizontal, as for the values up to n=16 it was already found in 
1951, and so we may again conclude that for all values of m higher than 
n=3 the ordinates are, roughly speaking, inversely proportional to n and, 
therefore, proportional to the horizontal dimensions of the topographic 
features of which the ordinates indicate the mean elevations. 

We shall now study the three curves of both figures in detail. In doing 
so we must realize that the greatest topographic feature is given by the 
continents and oceans and so both parts of the earth’s surface have this 
topography in common. As far as this great part of the topography is 
concerned, we may, therefore, expect that the sum of the squares of the 
topographic elevations, corresponding to a spherical harmonic term 
representing part of this great topographic feature, must be about the 
same in the land part as in the sea part of the earth’s surface, and, 
therefore, both equal to half of the sum of the squares over the whole 
earth’s surface. Assuming again that the land surface is 0.2917 times the 
whole earth’s surface we thus get the equation 


0.2917 12? = 0.7083 s? = 0.5 #2 
which gives 


(2) L=1.56s=1.31¢t 


We may expect that the ordinates of our curves, which are related to 
the continent-ocean boundary topography approximately obey these 
ratios. 

For the values of the three curves of both figures up to n=6 we can 
draw the same conclusions as in 1951. The first order term of fig. 1 is 
large, which corresponds to the well known fact that we can divide the 
earth’s surface into two hemispheres of which one contains a maximum 
of continent and the other a minimum, and that the difference is large. 
Roughly speaking we may state that the continents lay mainly on one 
side of our globe. One of the two centers of the big first order term is 
situated in the Black Sea south of Odessa and east of Constanza, and 
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this point is not far from the center of the land-hemisphere; the other is, 
of course, antipodic to it in the 8. Pacific. 

The ordinates of the 2nd—6th order terms form an upward wave, which 
shows a correlation to the wave formed by the Rayleigh numbers 
corresponding to convection-currents in the mantle of the earth distributed 
corresponding to the same orders of spherical harmonics. The smallest 
Rayleigh number corresponds to the greatest probability of the current 
distribution and so the inverse ratio of these numbers may serve as a 
measure for this probability. Table II lists the Rayleigh numbers An 
in the second column, the third column contains the ratio 600000/An, the 
fourth the value of ¢, the fifth the value of the third column multiplied 
by n#(n+1)', and the sixth the value of 7. The hatched curves in fig. 1 
and 2 show the values of columns 3 and 5 respectively. 


TABLE II 
ain An 600000/An t (600000/A,) x nt (mn + 1)? iM 
l 1722 | 348 979 493 1385 
2 825 727 841 1780 2060 
3 663-905 905 3135 | 3135 
4 683 879. | 875 3920 3979 
5 800 | 750 | 815 4110 | 4464 
6 1000 600 ~=~——-509 3877 3299 
i 1288 | 466 | 481 3487 3600 


Table IT and the curves in the figures 1 and 2 clearly show the correlation 
between the values of ¢ and 7 and the probability of mantle currents 
distributed according to the corresponding value of n. 

As the writer has previously discussed!), it is likely that we may 
explain these facts about the first six ordinates of both curves by the 
following hypothesis about the beginning of the earth’s history. Assuming 
that a period occurred when the earth was an undifferentiated fluid body, 
the cooling at the outside must have caused a first order convection 
system, which brought the heavy metals in the core and light sialic 
matter to the surface. The same current must have pushed the latter 
together in an urcontinent which had a surface about equal to the sum 
of the present continental surfaces. The presence of the core must have 
brought this current to a stop. When the cooling again caused convection, 
it adapted itself to the mantle and thus assumed a 3rd, 4th and 5th 
order spherical harmonic character. It tore apart the urcontinent which, 
being lighter than the mantle, could consolidate in the time-interval 
between the two convection-systems. The mantle matter, which by 


1) Venine Meresz, F. A., A remarkable feature of the earth’s topography, 
etc., Proc. Kon. Ned. Akad. v. Wetensch., Amsterdam, Ser. B. 14, 3 (1951); Venrnae 
Mezrnesz, F. A. and W. A. HersKaNEN. The Earth and its Gravity Field, McGraw- 
Hill Company, 1958, pp. 423 es., 11-6, 7, 8, 9. 
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consolidation got denser, subsided and melted again. So the mantle 
surface remained liquid during a much longer period and the parts of 
the urcontinent could, therefore, by the second convection current system 
be transferred to the areas above the subsiding columns of these currents 
and assume the present continent distribution. 

Examining the ordinates of the Ist to the 6th order, we find that the 
ratios of 1, s, and t, and, therefore also, of Z, S, and 7 check well with 
the values given by formulas (2), which is in harmony with the supposition 
that they indeed mainly concern the continent-ocean boundary topography. 

The further development from n=7 to n=31 shows a remarkable 
series of waves and maxima. Notable are the following values. 
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n= 7, peak only in ocean floor, 

n=10, peak about equal in size in continents and oceans, 

n=14, peak stronger in the oceans than in the continents, 

n=17, peak very strong in the continents and strong in the oceans, 
n=21, small peak only in the continents, 

n= 23, depression only in the oceans, 

n= 26, strong peak in the continents and smaller in the oceans, 

n= 28, slight peak, only in the oceans, 

n=30, strong peak in the continents and smaller in the oceans. 


Before discussing these results we have to examine whether there is 
evidence of multiples of order numbers which may, therefore, be connected 
with the same topography system. Multiples of n=3 and n=4 do not 
seem to be present. For n=5 the conclusion can not be quite so sure; 
the peak for n=10 may be connected with that for n=5, but there is 
not much evidence for higher multiples; it does not seem likely that the 
peak for n= 30 can be considered as such. Part of the peak for n= 14 may 
well be connected with that for n =7; the last occurs practically exclusively 
in the oceanic curve but for n=14 the continental curve also shows one, 
though decidedly smaller. For n=21 the peak in the oceanic curve is 
negligible and so it does not seem very probable that the small peak 
for n= 28, exclusively occurring in the oceanic curve, can be attributed 
to a multiple effect of the peak for n=7, though it is possible. The great 
peaks for n=17 and n=26 are evidently no effect of multiples and 
probably neither that for n=30. 

The peak in the oceanic curve for n=7 does not occur in Prey’s 
development of 1922 and so it must be caused by more recently discovered 
topography in the ocean-floors. Among these features a great part is 
probably played by the Mid-Atlantic and the Mid-Indian Ocean ridges. 
In the following way we can understand their having a 7th order character. 
We may imagine their cause to be connected with the great mantle 
convection system mentioned above, which broke apart the urcontinent 
and thus formed the Atlantic and Indian Oceans. These oceans must 
then have been located above the rising columns of this current system. 
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We can understand their main effect to have been a widening of the 
gaps in the urcontinent which they brought about. Still it seems explicable 
that in the middle part of the crust above the rising columns an amount 
of the sialic matter could remain present, which was not carried away 
by the currents towards both sides of the convection system. If this is 
indeed the origin of these mid-ocean ridges, the spherical harmonic 
development must represent them as terms of double the order of the 
terms corresponding to the distribution of the current system itself. 
As among these last terms those of the 3rd and 4th order played a large 
part, this would lead to these mid-ocean ridges having a 7th order 
character. 

In view of the dimensions of a few of the greatest island-ares, as e.g. 
the Indonesian are as a whole, including the Philippine archipelago, the 
Bonin, Mariana, Guam, Yap, Palau islands arc, and the Aleutian islands 
arc, it appears possible that they belong to the ocean floor topography 
represented by the 14th order term of our spherical harmonic develop- 
ment. 

The peaks for n=26 and n=30 may probably, at least partly, be 
attributed to the strong topographic effects of smaller island-ares, of 
which the sizes correspond to values of » in that neighbourhood. For a 
few of these ares, as e.g. for those inside the Indonesian Archipelago, in 
the Caribbean area and in the Mediterranean north of Corsica, there is 
strong evidence for the recent subsidence of the deep basins and the 
most likely explanation for these spectacular phenomena is given by the 
hypothesis of convection-currents in the upper layer of the mantle, which 
reach down in the transition-layer. The cause of such currents must 
again be looked for in the earth’s cooling at the surface, and the trigger 
effect needed for getting them into movement may probably be attributed 
to the concentration in the are around, because of crustal folding and 
downbuckling, of sialic rocks which are richer in radio-active constituents 
than the mantle rocks on which the sial floats. 

So here again we are brought to attribute the features revealed by the 
spherical harmonic development to mantle currents and so it seems 
indicated to explain also the other peaks, as e.g. the important one for 
n=17, to the same cause. 

These suppositions are supported by the following investigation, which 
shows that the peaks for n= 10, 14, 17, 21, 26 and 30 appear to correspond 
to currents in 1/2, 1/3, 1/4, 1/5, 1/6 and 1/7th parts of the thickness of 
the mantle. We shall first discuss the evidence in favour of this hypothesis 
and afterwards discuss it more closely. 

We can make two suppositions about the way the mantle divides 
itself in layers of current systems. We can assume the layers to have 
equal thickness, which we shall denote by de, or we may suppose the 
thickness of each layer to be proportional to its radius; in the latter case 
we shall indicate the thicknesses of the upper mantle layer by dp and for 
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the lowest one by d,. The values of dp, dy» and d, are listed in the 
columns 2, 5 and 6 of Table III. Column 3 of that table gives the half 
wave-length of the topography corresponding to the above mentioned 
values of n; for this purpose the writer has chosen the half wave-length 
L, on the equator for the sectorial sub-term, i.e. the value of 40.000/2n km; 
this value is, however, reduced to the middle plane of the current by 
substituting for 40.000 km the mean of this value and of the equator- 
length at the lower boundary of the current-layer. Column 4 gives the 
ratios re of these half wave-lengths L, divided by the thicknesses d, of 
the current-layer given by column 2, and column 8 the ratios ry of the 
half wave-lengths Ly for the equators, reduced to the geometric mean 
values between the surface and the lower boundary of the current-layers 
divided by the thicknesses dy; these mean half wave-lengths Ly are 
listed in column 7. 


TABLE III 

rr 
a de Le . | eS ies . 
km km | km | km = =km 2 
5 | 2900:1 = 2900 | 3090 1.07 2900 2900 2953 1.02 
10 | 2900:2 = 1450 | 1774 1.22 | 1668 | 1231 1720 1.03 
14 | 2900:3= 967 | 1320 1.36 | 1168 | 779 1292 iim 
17 | 2900:4= 725 | 1110 | 1.53 | 898 569 1092 1.22 
21 | 2900:5 = 580 911 1.57 729 | 448 s98 | 1.23 
26 | 2900:6 = 483 740 1,52) 9 6137)) | 370" 5) © 732) 1), 119 
30 2900:7 = 414 645 1.56 530 | 315 | 640 1.21 


Studying table III we find that, irrespective of which of the two 
divisions of the mantle-thickness is true, inclusive also of solutions between 
the two, the results show regular and acceptable values for the ratios r 
between the mean breadths L and the heights d of the convection- 
systems in the different layers. We must emphasize that it is not necessary 
to assume that the convection in the deeper layers is as active as in the 
top layer. This question probably depends on the way the temperature 
gradient changes from top to bottom of the mantle. 

The values of r, as well as those of rp show a curious increase between 
n=10 and n=17, while for the values n=5 and n= 10, as well as for the 
group of values for n=17, 21, 26 and 30, they are remarkably constant. 
This result is in good agreement with our hypothesis that all peaks of 
our spherical harmonic development curve have the same cause, viz. 
convection cells in the mantle. The increase of 7, and ry between n= 10 
and n=17 can probably be explained by the fact that for n=17, 21, 26 
and 30 the currents reach down to the transition-layer; this layer must 
facilitate the lower horizontal part of these currents, because part of the 
matter in the rising current is provided by the change, because of the 
current, of the mantle constituent, (Fe, Mg), SiO.s, from the heavier, 
cubic, in the lighter, orthorhombic, phase and the resulting increase of 
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volume, while part of the matter of the subsiding column disappears 
in the transition-layer because of the reverse phase-change. So the lower 
horizontal part of the convection-current is weaker than it normally 
would be and, therefore, meets less resistance. We can, therefore, under- 
stand that the horizontal dimensions are larger than they would otherwise 
be, and consequently the ratios r likewise. 

We may add that the values n= 14, 17, 21, 26 and 30 all correspond to 
the topographies of larger and smaller island-are basins and that the 
fairly rapid subsidence of these basins to great depth can best be explained 
by the relatively quick crystal phase-changes in the transition-layer. It is 
interesting to find that the dimensions of the basins appear to have a 
series of preferred sizes, given by the peaks in the spherical harmonic 
development curve. 

The basins are surrounded by geosynclinal belts under which the 
rising columns of the convection-currents are located; we supposed that 
the sial-concentration in these geosynclinal belts caused by the folding 
and downbuckling of the crust, which by its greater radio-activity brings 
about higher temperature, provides the trigger effect for the convection- 
currents. We may, therefore, also state that the peaks for n= 14, 17, 21, 
26 and 30 must correspond to the distribution over the earth’s surface of 
the present geosyncline belts. This leads to the supposition that the 
topographies of the great folded mountain ranges also contribute to 
these peaks. This checks with the fact that the ratio of the peaks in the 
oceanic and continental curves show greater figures for //s and L/S than 
those of formula (2) and especially so for the last three values of n 
mentioned, of which the peak for n = 21 is even nearly entirely continental. 
The continents, therefore, contribute more to these peaks than the oceans, 
and as geosynclines in the oceans can, eventually, only bring about 
low ranges of less than 1000 m high), our result corroborates the 
supposition that the high folded mountain ranges in the continental 
geosynclinal belts have considerably contributed to these peaks of our 
curves. 

A point which deserves to be discussed, is whether we can expect 
that for a small type convection-cell of a depth of e.g. 1/7th of the 
thickness of the mantle, or slightly more (for dp), this depth can be 
governed by the lower boundary of the mantle, which seems to imply 
that all the other six mantle layers must come into movement. For 
answering this question the writer thinks we have to realize that probably 
the whole mantle is crystalline and that we thus must assume the currents 
to have the character of permanent deformation of crystals. Two facts 
are in good harmony with this hypothesis. In the first place the velocity 
of mantle currents must be very small, viz. of the order of at maximum 
a few inches per year. Secondly the mantle matter appears to behave 


1) See e.g. HEISKANEN and VENING Mernezsz, The Earth and its Gravity Field, 
McGraw-Hill Company, 1958, p. 343. 
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plastically, i.e. for a small stress deviator it reacts elastically and the 
deviator has to exceed a certain elastic limit before permanent deformation, 
leading to flow, can occur. It is well known that this is a property, 
pertaining to crystals, which even for high temperature does not disappear, 
and this is in good harmony with the fact that even for the high temperat- 
ures at the bottom of the mantle we must assume this elastic limit to 
be present: if this were not the case, this inner part of the mantle would 
be continuously subject to convection, and this cannot be true as we 
must assume that in the same way as the whole mantle, it only does 
so during orogenic periods. 

This behaviour of the mantle can probably solve our problem. When 
the temperature gradients in the upper layer of the mantle are leading 
to a small type convection-cell, the flow is preceded by an elastic 
deformation and this deformation must have a vertical wave-length 
determined by the whole thickness of the mantle, divided by an integer, 
in case the mantle were a plane layer, and for the actual shape of a 
spherical shell, probably divided according to one of the types of division 
as given by our two cases, or according to one between the two. In any 
case the lower boundary of the mantle must play a part in this elasticity 
problem and so we can understand the discrete peaks in our spherical 
harmonic development; it is clear that the elastic deformation pattern 
must also be the pattern of the current system coming into being when 
the stress deviator increases above the elastic limit and leads to permanent 
deformation and to flow. In a subsequent paper the writer hopes to be 
able to give the solution for the problems here involved. 


9 Series B 


GEOPHYSICS 


THE RESULTS OF THE DEVELOPMENT OF THE EARTH'S 
TOPOGRAPHY IN SPHERICAL HARMONICS UP TO THE 
31st ORDER; PROVISIONAL CONCLUSIONS 
Il. THE PLASTIC-ELASTIC DEFORMATION AND FLOW IN THE MANTLE 


BY 


F. A. VENING MEINESZ 


(Communicated at the meeting of February 28, 1959) 


Before taking up this problem, the writer may point out that the type 
of smaller convection cell in the mantle, as derived from the spherical 
harmonic analysis of the earth’s topography, differs from what the author 
previously surmised and published in ‘“‘Gravity Expeditions at Sea’, 
Vol. IV, pp. 46 e.s.!) and, shorter, in ‘“The Earth and its Gravity Field”, 
pp. 4388 e.s.2). Basing himself on the supposition that the temperature 
gradient 6 diminishes downward as given by the formula f=foe~™ 
(z being the vertical coordinate, positive downwards), and assuming, for 
enabling him to deal theoretically with the problem, that the mantle 
consists of a Newtonian liquid of extremely high viscosity, he arrived 
at a solution according to which the small type convection-cell continues 
far down with a slow dying out downwards of the lower horizontal limb 
of the convection current. The difficult mathematical solution of the 
hydrodynamical equations was provided to him by Dr. J. G. VAN DER 
Corput (Amsterdam, now in Pasadena), who gave an elegant treatment 
of this complicated problem 8). 

The results now derived from the spherical harmonic analysis of the 
earth’s topography show a different picture, viz. a vertical distribution 
of many different convection cells in successive layers downwards in the 
mantle. For the 30th order peak in the development curve we even have 
to assume seven cells, one on top of the other, over the whole thickness 
of the mantle. The cause of the error in the picture formerly given, is 
easily traced in the fact that the mantle does not behave as a Newtonian 
liquid, but as a crystalline layer having a limit of elasticity and only a 
possibility of slow flow if the stress deviator surpasses this limit. 

It is difficult to give a complete quantitative treatment of this type 
of deformation, and so the writer shall limit himself to a geometrical 


1) Publication of the Netherlands Geodetic Commission, Delftsche Uitgevers 
Mij., Delft, 1948. 


2) McGraw-Hill, New York, 1958. 
3) Panta Rhei, “Simon Stevin’’, 1947. 
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analysis combined with a geophysical reasoning which seems acceptable. 
We shall start with the problem of a plane plastic layer subject to a 
small type of convection with a vertical division in several layers of 
cells, supposing adjacent cells—in a vertical as well as in a horizontal 
sense contiguous—to rotate in inverse sense. 

We shall treat the problem as a two-dimensional one with the coordinate 
x in the horizontal sense and z in the vertical one, positive downwards. 
We shall denote the elastic displacement in the x direction by w and in 
the z direction by w, and indicate the velocities in these directions by 
u and w. For wu and w we can put 


(1) u= Aq sin fx cos qz w= —Af cos fx sin gz. 
It follows 
ou ow 


So we see that the proposed elastic strain has a purely deviatoric character 
and that we, therefore, can use the same equations for the velocity comp- 
onents wu and wz in case we study the plastic flow. Figure 1 represents the 
flow curves for this case in an area slightly larger than one cell. The 
dimensions of such a cell, i.e. the half wave-lengths in the horizontal 
and the vertical sense are given by Lz,=a/f and L,=2/q, respectively. 
The equations for these flow curves, which likewise are those of the 
elastic displacements, can easily be derived. 


i eeaoaen A 


Fig. 1. Flow pattern in plane layer. 
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If dz and dz are the components of an infinitesimal part of one of 


these curves, we have 
dz _w__ fcos fx sin qz 
qsin fx cos qz° 


dx u 


Separating the variables x and z and bringing the right member to the 


left, we obtain 
qeosqzdz _ feosfadr _ 0 
sin qz sinfx 
and by integrating 


In sin gz+In sin fe=constant 


(3) sin fx sin qz=C 


which leads to the pattern of fig. 1. 

The equations (1) give the following results for the stresses, as well in 
the case of elastic deformations as for that of plastic flow. In the latter 
case we have of course to introduce % and w in stead of uw and w and yu 
then denotes the modulus of pseudo-viscosity instead of that of elastic 
shear. 


(4a) Oz = 2u i = 2u Afq cos fx cos qz 

ow 
(4b) Op? She 2u Afg cos fa cos qz 

OU : : 
(4c) Ty = ph (+3) = — uw A(q?—f?) sin fx sin gz. 


As for reasons of symmetry it may be expected, t, vanishes on the 
four sides of the cell (7=0, «=Lz, z=0, z=L,). Examining the formula (4) 
we may remark that for a certain magnitude of the deformation, 
characterized by the product fq, the deformation energy, which is 
proportional to the sum of the squares of the constituents of the stress 
tensor, here o%-+0)+2t?, becomes minimum when ty, vanishes, i.e. for 
f=q. If the phenomenon is free to adopt its cell distribution, we may, 
therefore, expect square cells. 


For our two-dimensional case, the equilibrium conditions are given by 


Oz , Ty 

(5a) St 4 St 4 Py = 
doz  oTy 

(5b) tae P,=0 


in which Pz and P; are the components of the force per unit of mass, 
which must work on the mass particle considered to give it the displacement 
given by the formula (1). By introducing the values for the stresses, 
given by (4) we obtain 


(6a) Pr= +nAg(g?-+f2) sin fe cos qz=p(g?+f2)u 
(6b) P,=—pAf(q?-+f2) cos fe sin qz=(g?-+f2)w 
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or, written vectorially 
(6) P=ule+ f2)V 


So the mass-force required to bring about the displacement as given by 
the formula (1) acts everywhere in the same direction as the displacement 
V and is proportional to it. 

We must realize, that if in part of the matter the limit of elasticity 
is surpassed, it will flow according to the same equations (1) in which, 
however, uw, w and wu are replaced by uw, w and 7, which denotes the 
modulus of pseudo-viscosity for this plastic flow. In this case the stresses 
will be the sum of the elastic and the pseudo-viscous stresses 
(tz = 2udu/dx + 2ndu/dx etc.). Even in case all the cells are undergoing 
plastic flow, the limit of elasticity can not be surpassed in the central 
parts of the cells, where the stresses are small or zero, and so these parts 
will only rotate and undergo a small elastic deformation. The deformation 
and flow pattern of fig. 1 will only be inappreciably modified. 

Before discussing these results, we shall take up the more difficult 
problem, presented to us by deformations and flow in the mantle of the 
earth, i.e. in a spherical shell. We follow a similar way as was hitherto 
followed and we shall limit ourselves to the study of the deformations 
and flow in the equatorplane while supposing a sectorial spherical 
harmonic distribution over the spherical shell. We shall further assume 
that the division of the mantle in layers is such, that their thicknesses dp 
are proportional to their radii in the same way as surmised in one of 
the two hypothesis of the preceding paper. We shall again suppose that 
the divergence /. V of the displacement vector is zero. We now have, 
instead of the components uw and w of this vector, the components v, 
in the east-west direction and v, in the direction of the radius @, i.e. 
in an opposite sense as w. We put 


(7) v,= —vA nsin nd V,=WoA cos nd 


in which v and wo are only functions of g, and A cos nd a sectorial 

spherical harmonic Y, of order nm in which, because we limit our 

deductions to the equator-plane, we have put the factor sin” #=1. In 

its general form the formula for v, equals (vo/sin J)d) Yn/dA and so, for our 
case, this becomes v9d¥ 7/2 1). 

For wo we put ae (- gk In $) 
wo = (n+ 1) 0? 

and, as the condition 7.V=0 provides us with the relation 


1 2) 
Oe ora) aaa i) 


1) See HEIskKANEN and Ventne Mervesz, “The Earth and its Gravity Field’, 


formulas (11-27), p. 413. 
2) Id. formulas (11-31), p. 414. 


we find 


We thus obtain 
‘ 4 1)A : 
(8) w%4= ue sin A cos (—aR In 2) Vo = — cos nA sin(—gR In 2) : 


These formulas have a similar shape as formulas (1) for the plane case. 
We can again replace v,, v, and mu by %,,%, and 7 for the case the limit 
of elasticity is surpassed and flow occurs. 

The values of 9 where v, vanishes, are given by 


(9) Os = Re ‘ak 


in which s is an integer: 1, 2,... (qR/z)In(R/r), where r is the inner 
boundary radius of the mantle. The maximum value sm gives the number 
of layers of convection cells in the mantle. From 


_7%),2 


44 it 


(10) i 


we can derive the value of gR corresponding to the number s, and, 
introducing s=1 in (9), the value of 01/R, which is the ratio of the inner 
and outer radii of all the cells corresponding to that number of layers; 
all cells are conformal. Putting R=6371.2 km and r=3471.2 km we find 


TABLE IV 
Sm qR ai/R (R—om)/dp 
1 5.173 0.54481 0.687 
2 10.346 0.73811 0.597 
3 15.519 0.81673 0.566 
4 20.692 0.85914 | 0.550 
5 25.865 0.88562 | 0.540 
6 31.038 0.90374 | 0.533 
7 36.211 0.91690 0.528 


The values of the third column, multiplied by R, lead to the thicknesses 
dy and d, mentioned in table III of the preceding paper. Fig. 2 gives 
a picture how, for a case of seven layers of convection cells and for 
sectorial spherical harmonic distribution of order n=30, the equator- 
plane is divided in cells. The horizontal dimension of each cell at the 
equator is given by aR/n. 


Fig. 2. 
Pattern of conformal convection-cells for i 
seven layers (for 1/10th part of the mantle). 
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For deriving the displacement-curves or flow-curves, respectively, we 
can follow the same line as for the plane case. The components of an 
infinitesimal part of these curves are odd and do and we get 
(n+ 1) cos (nA) sin G qR \n 2) 


cs es 
odA VA 


qR sin (nd) cos (- qk \n 2) 


Separating the variables, we obtain 


qR cos (- qR In 2) do 


(n+1) cos (nd) da 2 


sin (nA) 


sin (—qRIn 2) 9 


and, by integrating, 
(n+1) 


n 


In sin (nA) + In sin (- qh In 2) = constant 
or 


(11) en oe (- gR In £) =C. 


For one layer of cells, i.e. for sm=1 and for the order n=5 of the 
spherical harmonic distribution, i.e. for a breadth of the cells of 36° 
in the east-west direction, fig. 3 gives the pattern of the flow-lines. It is 
interesting to mention that this pattern checks with the curve for wo 
as function of 9, which the writer formerly derived for a general spherical 


Fig. 3. Flow pattern in the mantle for single convection-layer; equator-plane for 
sectorial spherical harmonic distribution, m = 5, 
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harmonic distribution of a convection-current over the mantle !); both 
deductions give a maximum value, which is situated considerably below 
the middle of the mantle cross-section, though according to our present 
deduction it is somewhat more excentric than according to the previous 
result. By differentiating v, with regard to 9 and equalling dv,/do to 
zero, we obtain the values of om where v,, resp. %,, has a maximum. 
We thus get 


(12) tg (—gR In) =—}qR 


from which om can be derived. For the ratio of the distance R— om from 
the maximum point M to the earth’s surface, to the height of the cell 
dp, we find the values listed in the fourth column of table IV. These 
ratios are valid for the corresponding distances and heights of all lower 
cells in the equator-plane which, as already mentioned, are all conformal. 

We may again point out, that for the case of flow, the limit of elasticity 
is not surpassed in the central parts of the cells, where the stresses are 
small or zero. So these parts will only rotate and undergo a small elastic 
deformation. The pattern of fig. 3 will only be slightly affected. 

As the deformations have a purely deviatoric character, the stresses 
causing them must be of the same type. We thus find 2) 


(13a) o,=2u (- “2 4 2) = 


0 oA 
6, (werd) 4 ngR 0 ; 0 
og Big tee te Sd t a) «gin ( — g 
2 ——s— 008 nd [4 cos ( gR |n 8) + sin ( qR \n 2) | 
3 On 
(13b) o.= 243? = — 


3 


Oy mina cos 22 [aR cos (- qk In 2) + 2 sin (- qR In 8) | 


(13d) y= (F242 2-2) — 


=u : sin nA [ (g2@R2— n(n + 1)) sin (- qR In 2) —3qR cos (- qR In 4) | 
(13e) T%=T,= 


For the case of flow we get the same equations but we have to replace 
4, v, and yw by %,,%, and 7. The total stresses are the sum of the stresses 
for the elastic deformation and for the flow. 

We see that the shear-stress t, again vanishes for the vertical cell 
boundaries given by nA=0 and nA=z, which, for reasons of symmetry, 
is clear. For the upper and lower cell boundaries given by —qR In 0/R=0 


1) The Earth and its Gravity Field, sec. 11-4, pp. 412 e.s. and Appendix to 
chapter 11, p. 442. 


*) The Earth and its Gravity Field, form. (11-42) in which # = 90°, p. 419. 
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and =z, this is only true for the first term of the factor between the 
square brackets; the second term does not vanish. The part of Tp 
corresponding to that second term equals — 3wv,/o. 

We can again study the minimum deformation energy, resp. minimum 
internal friction, for a certain magnitude of the deformation, now 
characterized by the product (n+1)¢R, and we find that for values of 
n=5 it is closely approached by the condition 


(14) q?@h?2=n(n+1) 


We may roughly state that the convection-cells may be expected to 
have a vertical dimension nearly equal to, though slightly less than, 
their mean east-west dimension. This is in good agreement with the 
results derived from the spherical harmonic development of the earth’s 
topography, as given in the last column of table II] of the preceding 
paper and as discussed in the ensuing paragraph. We find the same 
ratios if we divide the values of qR given in the first column of table IV 
by the numbers n of the corresponding peaks of the development curve 
(5, 10, 14, 17, 21, 26 and 30), and we obtain only slightly smaller values 
if we divide them by the square root of n(n+1). 

We may conclude that the good agreement is a valuable support for 
the mantle convection hypothesis developed here and for the adopted 
shapes of the convection-cells. If these cells were not conformally 
diminishing in size downwards, the upper cells, the horizontal dimensions 
of which are given by the spherical harmonic development, could not be 
nearly square in cross-section, as we have seen we must expect them to be. 

By introducing the stresses, as given by (13), in the equilibrium cond- 
itions, we can derive the components P, and P, of the mass-forces, which, 
for explaining the plastic deformation and convection-currents according 
to the formulas adopted here, must work on each particle. We obtain 

BOOE aig tt 2 OO 

g oA e@ %@ 
=u = sin nA [ oR {72k +n(n—1)} cos (- qgR \n $+ 


(15a) P,= 


+ {3 @2h2 + 2n(n + 1)} sin (—ak In al 


6, + 09— 206 009 1 Ot = 
PR) tis Tame gue 8 2g Tg Dh 


= Ih _ cos nA [ {(e+2) g2R2 + n2(n + 1)} sin (- qR |n 2) _ 
— (n+4) qf cos (—aR In 2)]. 


For large values of gR and n we may approximately write 


2R24 n2 
(16a) NV ear ee 


2R21 n2 
(16b) | nc ern ee 
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or 
(16) Pw pu LS V 
which corresponds to formula (6) for the plane case. So it remains 
approximately true that the required mass-forces have the same direction 
as, and are proportional to, the displacement vector. It is, however, 
important to notice that, as it follows from (15) as well as from (16) —if 
we realize that V is inversely proportional to o?—that P is inversely 
proportional to 04. This seems to indicate that the energy needed for a 
convection-current of the multiple layer type is in a greater measure 
due to the temperature gradient near the lower boundary of the mantle 
than to that near the earth’s crust. 

This result need not surprise us. We may assume that during the 
tertiary a great size convection-current has turned over the whole mantle 
and has brought the lower temperature matter down and into contact 
with the core, and the higher temperature matter up and into contact 
with the crust. In both contact zones a strong temperature gradient must 
have resulted, which we may expect to have led to convection of small 
size. The results found by the spherical harmonic analysis of the earth’s 
topography appear to indicate that these two small size convection 
systems formed systems between and thus united in one multiple system 
throughout the whole mantle of small convection-cells, of which the 
horizontal and vertical dimensions were determined by the system at 
the surface. 

It is true that the mass-forces as given by (15) and (16) can not have 
been present throughout the whole current-flow, but for the subsiding 
and rising columns, the thermal expansion effect on the density of the 
temperature differences between both columns, caused by the vertical 
temperature gradient, combined with the fact that the temperature 
conduction is so slow that it can be neglected, brings about mass-forces 
proportional to the current velocity which closely resemble the mass- 
forces given by (15) and (16). The upper cell, moreover, must be expected 
to have come into being because of a trigger effect formed by the 
horizontal temperature gradient between the higher temperature below 
the island-are, under which we suppose the rising column of the convection, 
and the lower temperature below the basin inside the arc, below which 
we must assume the subsiding column. So for this cell the upper horizontal 
limb of the current is subject to mass-forces which likewise resemble 
those of (15) and (16). 

It is difficult to say how many layers of small cells have been caused 
by the temperature-gradient in the upper layer of the mantle, as 
mentioned above, and how many by the lower temperature-gradient 
layer. It is, therefore, difficult to say whether during the uniting in one 
system over the whole height of the mantle, there were in the middle 
one or more layers where the deformation did not come above the elastic 
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limit, or whether they will all flow; the latter seems the most probable. 

From our discussion it follows that for two reasons we may assume, 
that these deep island-are basins of preferred size can only originate in 
the geological period following on that, during which the great half turn 
convection current reversed the whole mantle. In the first place the trigger 
effect required can be provided by the horizontal temperature gradient 
attendant on the heat produced by the sial concentration by folding in 
the geosynclinal belt, which is richer in radio-active constituents than 
the more basic rocks. In the second place, it is only in this period, that 
a great temperature gradient is present between core and lower mantle. 

The question arises whether these deep island-arc basins can also develop 
when only a trigger effect at the mantle surface is present. We may 
probably assume that in that case no deeper rows of cells, reaching down 
to the core, can form. It is difficult to make sure about this conclusion, 
but the writer is inclined to suppose that in that case no basin develop- 
ment is possible. He bases this surmise on the following facts and reasoning 
regarding the Indonesian island-arc. 

The evidence of the earthquake foci at depths of 400-700 km below 
the deep basins checks well with the great values formerly derived !) of 
the shear-stress in planes under 45° with the vertical in the lower part 
of the subsiding column of the upper cell. This is in harmony with the 
supposition, that the deep basins of the archipelago are brought about 
by the small convection-cells already mentioned. The deep earthquake 
centers, however, continue much further west than the deep Banda basin ; 
we find them e.g. in the Java Sea north of West Java, although this sea 
is shallow. 

In view of our present discussions there seems to be reason for assuming 
that these shear-stresses, causing the deep earthquakes north of Java, 
are brought about by only one convection-cell layer, perhaps even moving 
with still lower velocity than the convection-currents here generally 
discussed. It appears possible that, because of the absence of a strong 
temperature gradient between core and lower mantle, no deeper cells 
did develop. The absence of the lower temperature gradient might be 
explained by the supposition that the great size mantle convection-cell, 
rising under Asia, has not yet proceeded far enough in its half turn to 
bring the low temperature mantle matter back below this part of the 
East Asian continent. If this surmise is right, the front of this lower 
flow back towards Asia is defined by the west and northwest border of 
the deep basins. The fact that mean elevation of the topography is 
here deep below sea-level, as well as the positive gravity anomalies in 
these basins, may probably be attributed to this great size current of 


1) See e.g. “The Earth and its Gravity Field”, 11-10, p. 438 es., where Fig. 
11-6 has to be corrected for the present supposition of cells, which are about square, 
but for that case the strong shear-stress near C remains true. We can also derive 
this result from the stresses oz and oz of formula’s (4a) and (4b) of this paper. 
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relativily low temperature and, therefore, higher density. At the same 
time this current would cause the lower temperature gradient zone 
and thus bring about the series of convection-cell layers, reaching 
down over the entire height of the mantle. 


We shall, finally, take up our problem for a general spherical harmonic 
distribution of the convection system over the mantle. In this way we 
shall come to a still nearer approach to our problem of the origin of the 
different topographic features at the earth’s surface. 

For currents in a spherical shell we can put ) 


oYn , _ vo oY n ——— r 
ems) eee ae 9 a aaa ay} yl v.= wy, 


in which Y, is now a spherical harmonic of general character, and in 
which wo and vp are again only functions of @ which for 7. V =Oare related 
by the equation, already mentioned on page 4. Introducing the same 
formulas for wo and vo as we did there, we obtain 


R Oo on 
(18a) Vg = — £5 cos (—gR ng) 
_ qR o\ 1-2, 
(18b) Y= Pi (—qR In 3) ain 
> a (n+1) . Q r 
(18e) 0 = — a Bin (— aR In 2) bee 


For these more general formulas we can draw the same conclusions 
for the different values of gR and of 0;/R for different values of sm as 
given by table IV. We find the same equation as given in (12) for the 
radius for which v, has a maximum value and we thus find the same 
ratio’s of R—opm to the thickness dp of the cell-layer, as given by the 
last column of table IV. We find again that all cells, vertically below 
each other, are conformal; the spherical harmonic distribution of all 
cell-layers for a certain system of convection is identical. However, as 
long as we do not further define this distribution, we can not give any 
pattern of flow. That the conditions for minimum deformation energy for 
the elastic deformation and for minimum internal friction for plastic 
flow are again given by (14) can easiest be derived by making use of the 
fact that the formulas for the respective horizontal distributions of the 
current systems in the plane layer and in the spherical shell merge into 
each other, if we put 2) 


and as for the plane layer the minimum energy, resp. friction condition 
is given by q=f, this leads immediately to the condition (14). 
The conclusions we came to for the previous case and for the equator- 


plane, remain the same for the general case here dealt with and for the 
whole mantle. 

1) The Earth and its Gravity Field, pp. 413 and 414. 

2) Id. formula (11-32), p. 414. 


BIOCHEMISTRY 


STUDIES ON FRUCTOKINASE 


1. SUBSTRATE SPECIFICITY 
BY 


Cu. M. A. KUYPER 


(Communicated by Prof. H. G. K. WESTENBRINK at the meeting of Nov. 29, 1958) 


INTRODUCTION 


As substrates for fructokinase only the ketohexoses D-fructose, L- 
sorbose and D-tagatose have been described [1, 2, 3]. In nature several 
other ketoses have been found [4], and little is known about their meta- 
bolism in animals. Therefore it appeared interesting to investigate whether 
fructokinase can phosphorylate more ketoses and may be to study the 
specific stereochemical requirements for substrates. 


MATERIAL AND METHODS 


The enzyme was prepared from rat livers by homogenization in six 
volumes of 0.9 94 NaCl solution. The homogenate (I) was brought to 
Pu 4.5 with 3N HCl and then immediately to py 5.5 with IN NaOH [2]. 
After spinning down the agglutinated particulate fractions the clear super- 
natant (II) was heated for 10 minutes in a waterbath at 70°C. The 
precipitate was filtered off and the clear filtrate (III) fractionated in the 
cold with acetone. The fraction precipitating between 25 and 55 °% acetone 
was dissolved in destilled water and freed from acetone under reduced 
pressure (IV). These four stages are referred to in the following as Fractions 
Pro Le 

The last fraction contains 90 to 100 % of the initial activity of the 
homogenate, and is about 15 times as active per mg protein. Appropriate 
tests showed that it was free from interfering enzymes as adenosine- 
triphosphatase, hexokinase, fructose-1-phosphatase, fructose-1-phosphate 
kinase, aldolase and myokinase [5]. 

The enzyme assays were carried out in the following mixture: 1 mg 
fructose (or other sugar), 1.5 millimoles potassium acetate, 10 micromoles 
adenosine triphosphate (ATP) and 10 micromoles MgSO,; total volume 
1.5 ml. The amount of enzyme added is such that maximally about 
400 micrograms sugar are phosphorylated. Under these conditions the 
enzyme is saturated with the carbohydrate substrate, and its activity 
appears to be proportional to its concentration. The mixture was incubated 
during 30 minutes at 30° C, thereafter deproteinized by adding 2 mis. 
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ZnSO, and 2 ml 3N Ba(OH), [6]. After shaking the mixture was made 
up to 10 ml and filtered. From the filtrate an aliquot part was used for 
the determination of residual free sugar. Ketosugars with 5 to 7 C atoms 
were determined colorimetrically by the method of Ror [7], in the modi- 
fication of Hurs et al. [8], or titrimetrically by virtue of their reducing 
properties according to HAGEDORN and JENSEN [9]. The enzyme activity 
is expressed as micrograms free sugar disappearing in 30 minutes 
(= Aug sugar). 


RESULTS 1) 


1. Stereochemical specificity of the enzyme 

We first attempted to determine affinity constants Km for several 
sugars, but no linear relationships were found between the reciprocal 
values of reaction velocities and substrate concentrations (see LINEWEAVER 
and Burk [10]). Therefore only qualitative data are recorded in Table I ?). 
At the same time the interference of several sugars with the phosphoryla- 
tion of fructose was examined. Also a number of aldoses, disaccharides 
and sugar derivatives were tested as to their suitability as substrates or 
inhibitors: D-glucose, D-mannose, D-galactose, D- and L-arabinose, D- 
and L-ribose, D-xylose, D-fructose-1-phosphate, D-fructose-6-phosphate, 
x and f§ diacetonfructose, monoacetonfructose and sucrose. These com- 
pounds were not phosphorylated and did not display inhibitory action. 

It can be seen from Table I that only ketoses having a certain specific 
configuration at the first three C atoms of the sugars are phosphorylated. 
In general the configuration of the residual part of the substrate molecule 
appears to be immaterial; only bulky substituents (e.g. phosphate on C 6) 
abolish the possibility to react. The sugars that can be phosphorylated 
interfere also with the utilization of fructose, so there seems to be a common 
active site for all these substrates on the enzyme. 

According to Cort [11] only the furanoid form of fructose is phosphor- 
ylated by yeast hexokinase. We also examined the influence of the ring 
system on the rate of phosphorylation by working at 18° C, at which 
temperature the velocity of interconversion of the two forms is very low 
[12]. The determinations were performed in two systems: in the first the 
fructose was added in solid form at zero time, so that only the pyranoid 
form was present; the other contained an equilibrated solution of fructose, 
consisting of approximately 70 °% pyranose and 30% furanose [12]. 
After 15 minutes the respective amounts of fructose phosphorylated 
appeared to be 242 and 245 yg. If the pyranoid form were phosphorylated 
faster than the furanoid or vice versa these figures would have been dif- 


‘) These results form part of a thesis published in Dutch [5]. 
*) Iwish to express my gratitude to Prof. Reichstein, Basel, for gifts of D-tagatose 
(prep. B) and L-arabinulose, to Prof. Havinga, Leiden, for gifts of L-sorbose and 


D-tagatose (prep. A) and to Dr. Richtmeyer, Bethesda, for gifts of L-sorbose and 
ketoheptuloses. 
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ferent, so it seems safe to conclude that the form of the ring in fructose 
is of no importance for the action of fructokinase. 

The only exception to the stereochemical configuration requirements 
seems to be D-mannoheptulose. Although the configuration is adequate 
it is not phosphorylated, nor does it interfere with fructose utilization. 
It was considered possible that in solution this sugar does not possess a 
ring structure, but exists in an open-chain form. In accordance with this 
hypothesis the literature only mentions one value for the specific rotation 
of this sugar [4] (no mutarotation). We could further confirm this pre- 
sumption by determining various ultraviolet absorption spectra. Figure 1 
shows that D-mannoheptulose has an absorption maximum at 2650 A, 
indicating a free ketogroup [13], while D-fructose and D-galaheptulose, 
which possess no free ketogroup, have no absorption band with a maximum 
at this wavelength. 


TABLE I 
Substrate specificity of fructokinase 


The reaction mixtures contained 1 mg of the compound tabulated. ‘‘Utilization’’: 
D-fructose absent; “Inhibition”: 1 mg fructose also present; inhibition of phos- 
phorylation of the latter. n.t. = not tested 


Substrate Structure Utilization Inhibition 
: =O), 6 0 SN 
IDR) 2 5 6 a 6 6 & Iq Aes 5 << . ae = 
line A Ds 
t=) 
L-sorbose . b Ne SF sF 
A and B i -b + 
D-tagatose (prep. A an ) fs DA 
a 
IDWOpayNeeO) 2 5 5 5 5 CH,OHCCH,OH *) = 
0) 


L-arabinulose ..... . ek is n.t. 
7— 
D-mannoheptulose. .. . je a a nal 


0. 
D-glucoheptulose. . . . . Pie aie az 


ie 

L-galaheptulose .... . K 1% alle iastie 
0 

D-altroheptulose. . .. . iy: oa Taf 


*) According to LeurHarpT and Testa [3]. 
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Fig. 1. Absorption spectra in the near ultraviolet of D-mannoheptulose (A), 
L-galaheptulose (B) and D-fructose (C). 


2. The homogeneity of the enzyme preparation 

Of course it would be possible that our enzyme preparation contained 
a mixture of enzymes, each specific for only one or very few substrates. 
We examined this possibility by testing the relative activities of our 
preparations with D-fructose, L-sorbose and D-tagatose !) as substrates. 

In the various stages of purification of the enzyme the activity of the 
separate fractions was tested with the three sugars. The results (Table IT) 
demonstrate clearly that in the course of fractionation the ratio of the 
activities changes in respect to fructose and sorbose, while the ratio for 
fructose and tagatose is constant. Similar results of the same kind are 
shown in Tables III and IV, which give the activities in respect to the 
three ketoses of several samples of Fraction IV after storage at 2° C or 
heating respectively. Under both conditions the activity of the enzyme 
declines when assayed with sorbose, while it remains constant when tested 
with fructose or tagatose. This seems to indicate that sorbose on the one 


‘) The preparations of D-tagatose were obtained from two different laboratories. 
These samples, designated as A and B, are somewhat different in properties, and 
therefore also gives different results relative to fructose. 
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hand and fructose and tagatose on the other hand are phosphorylated 
by different enzymes. 
There are, however, also phenomena indicating that the same enzymatic 
entity is responsible for the phosphorylation of both fructose and sorbose: 
1. Intraperitoneal injection of either fructose or sorbose into rats 
causes an inductive increase of the kinase activity of the liver towards 
both fructose and sorbose (see Table V). In these experiments quantitative 


TABLE II 
Activity of Fractions I to IV with substrates D-fructose, L-sorbose and D-tagatose 
The activity is expressed as wg sugar phosphorylated by 0.5 ml enzyme solution, 
with the exception of Fraction IV of which only 0.3 ml was added 


Teton Activity Ratio of activities 
D-fructose L-sorbose D-tagatose *) fruct./sorb. fruct./tagat. 

I 289 209 116 1.38 2.49 

II 437 288 170 1.51 2.57 

TEE 461 268 188 1.72 2.46 

Ves) 473 303 186 1.56 2.54 

DVvae) 329 132 131 2.49 2.52 


*) Preparation A. 
**) Different preparations. 


TABLE III 
Activity of fructokinase preparations (Fraction IV) with substrates D-fructose, 
L-sorbose and D-tagatose *) after storage at + 2°C 
The activity is expressed as wg ketose phosphorylated by 0.3 ml enzyme solution. 
F = D-fructose, S8 = L-sorbose, T = D-tagatose 


Days Preparation I Preparation II Preparation III Preparation IV 
Bite eee LenS eee) IBY ys) GR TEs) IBY Mts) AS) mS) IRS 


0 403 198 2.03 386 244 164 1.58 2.37 382 192 1.98 473 303 1.56 
1 404 197 2.05 —- — - 

3 412 160 2.67 

5 404 114 3.54 

7 —_—- — — 398 164 162 2.44 2.47 380 0 ~ 475 0 ~ 


*) Preparation A. 


TABLE IV 


Activity of fructokinase with substrates D-fructose, L-sorbose and D-tagatose *) after 
heating 
The activity is expressed as wg ketose phosphorylated by 0.3 ml enzyme solution. 
Symbols for sugars as in Table III 


Minutes Activity Ratio of activities 
at 70°C F S) iD F/S E/T 
Ee SE eee 
5 Bye) ilalsy Uy 2.99 1.97 
10 325 83 165 Syl ia} 
15 335 Wes IMG 4.60 1.99 


*) Preparation B. 
10 Series B 
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extraction of the enzyme from the liver was effected by repeated homo- 
genization and centrifugation. The activity of the extract was expressed 
in ug ketose phosphorylated per 10’ nuclei in the homogenate [5]. 

2. The affinity for the magnesium-ATP complex [14, 15] is the same 
when Fraction IV is tested with either fructose or sorbose (see Table V1). 
The Km is 0.10 ml of a 0.10 M solution of this complex with either sugar. 

3. Both sugars compete for the same site on the enzyme surface (see 
Table VII). The amount of each sugar in the mixture was determined by 
combining the colorimetric and titrimetric procedures; the ratio of the 
molar extinction coefficients of the two sugars is different from the ratio 
of their reducing values. 


TABLE V 


Induction of the kinase activity towards D-fructose, L-sorbose and D-tagatose *) by 
D-fructose and L-sorbose 
The induction was effected by daily intraperitoneal injection during 7 days of | ml 
of a 40 % fructose or sorbose solution. The activity is expressed as ug ketose 
phosphorylated per 107 nuclei in liver homogenate. 
Symbols for sugars as in Table III 


Activity Ratio of activities 
Inductor Fr S 7. F/S B/T 
D-fructose 376 131 190 2.88 1.98 
L-sorbose 344 117 170 2.94 2.02 
Initial activity 220 76 111 2.90 1.99 
*) Preparation B. 
TABLE VI 


Affinity towards the ATP-Mg complex with D-fructose and L-sorbose as substrates 
The activity is expressed as ug ketose phosphorylated by 0.5 ml of enzyme solution 


a 


ml of 0.05 M ATP-—0.05 M Mg Activity 
solution D-fructose L-sorbose 
Se a 
0.10 195 — 
0.20 255 163 
0.30 296 190 
0.40 303 201 
0.50 — 210 
TABLE VII 


Activity of fructokinase (Fraction IV) towards D-fructose and L-sorbose in a miature 
of both sugars 
The reaction mixture contained 1 mg of each sugar. Activities are expressed as ug 
ketose phosphorylated by 0.3 ml enzyme solution 


. Activity 
Substrate D-fructose L-sorbose 
Exp. 1 Exp. 2 Exp. 1 Exp. 2 
D-fructose 598 520 — — 
L-sorbose — — 409 275 


D-fructose + L-sorbose 565 505 200 70 
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DISCUSSION 


From the results given above it can be concluded that “fructokinase”’ 
or “ketohexokinase”’ [3] might better be called ketokinase as not only 
fructose or ketohexoses but quite a number of ketoses are phosphorylated. 
From the rather extensive list of compounds tested it can be seen that the 
most important structural requirements for a substrate seem to be the 
following: 

1, The sugar must be a ketose. 

2. The sugar must possess a ring structure. Open chain ketoses are 
not phosphorylated. The type of the ring (pyranoid or furanoid) is 
immaterial, for both forms of fructose are phosphorylated at the same 
rate, just as the furanoid L-arabinulose [16] and the pyranoid L-gala- 
heptulose. 

3. The hydroxyl group at the third C-atom must be free and in trans 
position to the hydroxyl group at C 2. The fact that D-xylulose, which 
fulfills these requirements, can be metabolized by the rat, in contrast to 
the L form of this sugar [17], is in accordance with this conclusion. 

A very interesting feature is the incongruous behaviour of the enzyme 
towards fructose and sorbose. The active principles are certainly not the 
same, as in the course of preparation of the enzyme the activity remains 
constant when fructose is the substrate, while it decreases when tested 
with sorbose. Also the very variable ratios of activities found with both 
sugars suggest the occurrence of different active principles in the prepara- 
tion. On the other hand indications have been obtained suggesting that 
only one enzyme is responsible for the phosphorylation of both substrates: 
1. The same affinity towards the ATP-Mg complex, when determined either 
in the presence of fructose or sorbose; 2. The induction of “fructose 
activity’ not only by injections of fructose but also by injections of 
sorbose; 3. The induction of ‘‘sorbose activity” by injection of either 
fructose or sorbose; 4. The mutual inhibition of phosphorylation by both 
sugars, As a possible explanation of these conflicting data the concept is 
tentatively advanced of one enzyme molecule bearing two different but 
interacting active sites. One active site is concerned with the activation 
of fructose (and tagatose), the second with the activation of sorbose. These 
sites must be very close to each other, possibly they even have some 
structure in common, because fructose inhibits the reaction with sorbose 
and vice versa. Both sites cooperate with the site of attachment of ATP. 
The two configurations on the protein surface, responsible for the two 
activation sites do not have the same stability, since the rate of heat 
inactivation depends upon the choice of the sugar, i.e. whether fructose 
or sorbose is used for the determination of the degree of inactivation. 
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SUMMARY 


1. An enzyme preparation is made from rat liver, possessing activity 
towards fructose, that is free from interfering enzymes. The conditions 
of assay are such that a linear relationship exists between amount of 
enzyme and activity. 

2. The ability of this preparation to phosphorylate a number of 
sugars and sugar derivatives was tested. 

3. Of twenty three compounds tested, only ketoses were phosphoryl- 
ated. Common features of these ketoses are: a ring structure (pyranoid 
or furanoid) and a free hydroxyl group at the third C-atom in trans 
position with respect to the hydroxygroup at the second C-atom. The 
number of carbon atoms of the ketoses may vary from 56 to 7. 

4. The name ketokinase is proposed for this enzyme instead of the 
older names fructokinase or ketohexokinase. 

5, Arguments are given for the tentative conclusion that the activity 
towards the various sugars is caused by only one enzyme. This enzyme 
seems to possess at least two different active sites for the several sugars. 


From the Laboratory of Cytochemistry, 
University of Nijmegen 
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ERRATUM 


In the paper by M. P. H. WEENINK and P. Groen, “‘A semi-theoretical, 
semi-empirical approach to the problem of finding wind effects on water 
levels in a shallow partly-enclosed sea’’, Proceedings Koninkl. Nederl. 
Akademie van Wetenschappen, Ser. B, 61 (No. 3), 1958, pp. 198-213, 
on page 206, formula (24) is not correct; it should read: 


( A(HH) = {(1.03 — 0.75 x 10-8 A) cos « — 
i ee — (0.10 + 0.26 x 10-6 2) sin «} 9-1g-1 Dy a f(x) Vi, 
where a= width of the North Sea model =5.75 x 105 m and D;= charac- 
teristic depth of area [=43 m. 


As 


ect 
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